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Abstract 

We derive some basic equations related to the redshift drift and we show 
how some dark energy (DE) properties can be retrieved from it. We consider 
in particular three kinds of DE models which exhibit a characteristic signa- 
ture in their redshift drift while no such signature would be present in their 
luminosity-distances: a sudden change of the equation of state parameter wde 
at low redshifts, oscillating DE and finally an equation of state with spikes at 
low redshifts. Accurate redshift drift measurements would provide interesting 
complementary probes for some of these models and for models with varying 
gravitational coupling. While the redshift drift would efficiently constrain mod- 
els with a spike at z ~ 1, the signature of the redshift drift for models with 
large variations at very low redshifts z < 0.1 would be unobservable, allowing 
a large arbitrariness in the present expansion of the universe. 
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1 Introduction 



Data suggests that our universe has entered a stage of accelerated expansion rate. 
This is a radical departure from conventional cosmology in which the expansion was 
constantly decelerated except for the inflationary stage in the very early universe. The 
mechanism behind this late-time accelerated expansion is unclear and many models 
were suggested [T]. If it is caused by some isotropic perfect fluid component called 
dark energy (DE) while gravity is governed by General Relativity (GR), then DE 
should make up about two thirds of the universe content. The simplest possibility 
a cosmological constant A is problematic from the theoretical point of view because 
of its tiny amplitude. Further while the data is in agreement with a cosmological 
constant A, the precision of these data do not allow to rule out models where the 
equation of state (EoS) parameter of the dark energy component wde differs from — 1 
and varies with time, and in which the corresponding energy density depends on time 
as well. A great number of models have been proposed and all their properties were 
intensively investigated in the hope that future accurate observations would allow to 
discriminate between these models and that only a small part of them would remain 
as viable candidates while the other models would be ruled out. Another attractive 
possibility is a modification of gravity on cosmic scales. Virtually all the foundations 
of standard cosmology have been questioned in the quest for a solution to the dark 
energy problem. Decisive progress in settling the issue will come from observations. 
Experiments of different kinds are planned that will probe with exquisite precision 
the universe background expansion and also the evolution of the perturbations. 

To make progress in the understanding of the nature of DE it is desirable to 
explore all ways in which its properties can be probed. One such probe is the redshift 
drift. The feasibility of such measurements was reassessed in [2], many years after 
its theoretical discovery [3]. It has been suggested recently as a way to test general 
properties of the DE paradigm |4j. These measurements could give additional insight 
into the properties of DE because they probe directly the quantity H(z). It is this 
aspect that we want to emphasize and to investigate in this work. Before embarking 
on the quantitative assessment of redshift drift data for some specific models we first 
derive some general properties of the redshift drift relevant for the study of DE. 



In an expanding universe many physical quantities evolve with time. The redshift 
suffered by radiation emitted by a source depends on time as well. Indeed let us 
consider radiation emitted by a source at the emission time t e . If this radiation is 
observed at time to, the corresponding redshift is given by 



where we use the notation z(t e ,t ) to emphasize that the redshift depends on the 
emission time t e and on the observation time to- If radiation emitted by the same 



2 The redshift drift 




(i) 
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source is observed at time to + Sto, the redshift will change by an amount Sz 

Sz = z(t e + St e , t + St ) - z(t e , t ) . 
It is easy to derive the following expression 

Sz = — (a - a e ) 
a e 



H 5t {l + z{t e ,t )-h e ) 
H St (l + z- h{z)) 



he = JT 
-"0 
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(4) 
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with the obvious notation H e = H(t e ). We have dropped the subscript e and finally 
we return to more conventional notations setting z = z(t e ,t ). It is obvious from ([3]) 
that 5z is a decreasing negative function of 2 in a universe whose expansion rate is 
always decelerated. This is what would happen e.g. for an Einstein-de Sitter universe. 
However the situation changes when at least part of the expansion is accelerated. Let 
us consider a universe for which the deceleration parameter q satisfies q < for z < z q 
and q > for z > z q and q(z q ) = 0. This is the case for a flat universe with dust-like 
matter and a cosmological constant A satisfying f^o > §• For this universe we see 
from ([3]) that 5z must be an increasing function of z on the interval z < z q during 
which the expansion rate is accelerating. At large redshifts, when the expansion rate 
is decelerating (and matter-dominated), 5z will decrease in function of z. Note that 
Sz vanishes in a Milne (empty) universe. 

It is interesting to study this behaviour by considering the slope (Sz)' = 4- (Sz). 
We have 

(Sz)' = H 5to (1 - ti) . (6) 
The slope is positive for hi < 1 and negative for hi > 1. For a universe with accelerated 



expansion rate on the interval < z < 
h'(zm) = 1. One can show that we have z m 



Sz reaches its maximum at 



> z q with z r , 



when 



typically only slightly larger 



than z q . We will consider below even more sophisticated scenarios where the slope 
is first negative at z = 0, then positive at higher z in order to ensure a late-time 
accelerated stage and then again negative, 
dimensionless quantity 



It is convenient to introduce the new 
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It is easy to derive the following equalities 
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where spatial flatness is assumed in ([9]). When we use luminosity distances, second 
derivatives are necessary. We have in particular 



-qo 



(10) 
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In contrast to the luminosity distance (Il(z), the leading order of Sz for z <C 1 is model 
dependent and allows to discriminate between models with same H Q but different go- 
Physically this is so because the slope at z = is sensitive to the dark energy EoS 
as we see from eqs. (15f§|) . Note that at the redshift z\ ^ where 5z = (if it exists), 
(Az)' = —q. The following inequality 

(A*y > i , (ii) 

is satisfied by a universe on all redshifts for which it is in a phantom phase, q < — 1 or 
Weg = VL DE wde < — 1- Hence a universe in a phantom phase on small redshifts must 
satisfy inequality (fTTj) on these redshifts. If it is DE itself which is of the phantom 
type, wde < — 1, the following weaker inequality has to be satisfied 

{Az)' > 1 - l ^^+f . (12) 

When DE is a cosmological constant A, the above inequality becomes a strict equality. 
Note that (Az)' = 1 in a de Sitter universe which corresponds to the asymptotic future 
of a flat universe with a cosmological constant A (Q m — > 0). 
These properties can be seen with the small- z expansion 



Az = 

which is to be contrasted with the expansion 



(13) 



D L (z) = H d L (z) = z + i(l - q )z 2 + ... (14) 

The leading term of ffl4|) is the same for all models with identical Hq in contrast 
to ( TL3]) . The result (|T0|) is obviously recovered from (TlBl . In (fT3|) j is the "jerk" 
parameter, j = at z=0. Expression ( !T3|) shows also how Az deviates from a 
linear law at small redshifts. From ( fT5|) . this deviation is positive i.e. 

Az > -q z , z > 0, z w (15) 

on some small interval around 2 = when the acceleration is slowing down on small 
redshifts. These properties of the redshift drift are illustrated with Figure 1. 

Let us return finally to the amplitude of the redshift drift. Clearly it is proportional 
to the small quantity H 5t . In other words all equations and inequalities derived 
above make use of Az, the redshift drift in units H 6t . For 6t = N years, H 6t = 
1.023 N hx 1(T 10 (here h = H / 100km/ sMpc), so 5t should be at least of the order 
of 10 years to yield a measurable effect. We will consider in the next section some 
specific models. 



3 Models with peculiar equation of state 

We consider now various DE models exhibiting a characteristic signature in their red- 
shift drift because the equation of state (EoS) parameter wde has a special behaviour 
on low redshifts. The choice of models in this study looked for distinctive features. 
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Figure 1: a) On the left panel, three models are displayed with a large variation 
of wde on low redshifts. In these universes the acceleration is slowing down on low 
redshifts and the expansion is even decelerating on very low redshifts for the top (pink) 
and bottom (green) curves. The red curve (in the middle) is taken from [8] and has 
go ~ 0. b) On the right panel the corresponding redshift drifts are displayed. The 
blue curve corresponds to a ACDM universe. Note that the other three curves satisfy 
the inequality ( |T5l) . The corresponding luminosity distances have no characteristic 
signatures though the pink curve is in tension with SNIa data while the model with 
the green curve provides an excellent fit because strong variation of wde for this 
model, with presently decelerated expansion, is confined to the interval z < 0.1. 

3.1 Large variation of wde on low redshifts 

We start with models with a large variation of Wde starting from very low redshifts 
and on until today. There are several situations where such a variation can occur, 
for example in models where the accelerated stage is transient and the accelerated 
expansion is already slowing down on very low redshifts. This can be the case for 
some quintessence models (see e.g. [5], [6]) or for more exotic models with a singu- 
larity in the future (see e.g. [?]). It was also proposed recently as an interesting 
phenomenological ansatz that can improve the fit to SNIa data compared to a pure 
cosmological constant [S|. Three models are shown on Figured! We see that these 
models have a characteristic slope of their redfshift drift at very low redshifts which 
differs significantly from the slope obtained for ACDM (Az)' = 1 — f ^+i-Az* • ^ ne 
slope at z = is negative if the universe is presently decelerating and rapid variation 
of Wde produces a sharp change in 5z. Redshift drift measurements on such very 
low redshifts suffer from increased difficulties related to the peculiar velocities of ob- 
jects and are virtually impossible. While the signatures of the models considered in 
this subsection are most probably not observable, they illustrate nicely the results 
derived in section 2. As a substantial amount of SNIa data have been added in the 
interval 0.1 < z < 0.4 in the Union2 set as compared to the Constitution set, two of 
these models are now in tension with SNIa data. The model with the pink curve on 
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Figure 2: a) The redshift drift for several models with a "spike" in their EoS pa- 
rameter wde is shown on the left panel. The spikes located at z — 0.2 (red and 
green curves) have identical widths and amplitudes (maximal heighths) A sp = 0.5 
and A sp = 0.2 respectively. The third spike at z = 1 (pink curve) has a larger width 
and amplitude A sp = 0.7. The location and amplitude of the spike leaves a neat 
signature in the redshift drift, b) The redshift drift of ACDM and a model with 
oscillating wde (with amplitude A osc = 0.3) are displayed on the right panel. The 
oscillations are clearly seen in the redshift drift while they are essentially erased in 
the luminosity distance. 

Figure 1 (denoted mod2 in Table 1) is in tension with the data because it departs 
from Woe — — 1 on higher redshifts z < 0.5. However the model with the green 
curve in Figure 1 (denoted modi in Table 1) provides an excellent fit to the data. It 
has a sharp departure from woe — ~ 1 on redshifts z < 0.1 only and a decelerated 
expansion today. This is in agreement with results obtained in [5]. 

3.2 Spikes in the equation of state 

Here we consider models where wde exhibits a spiky feature (see e.g. [9], [10]). The 
larger the redshift where a spike of given amplitude (i.e. maximal heighth) A sp , and 
width is located, the smaller the effect on 5z because Qde decreases with redshift and 
hence the effect on H(z) of a spike in wde will decrease as well. We consider models 
which have the same fiducial value wde,o = —1 with a superimposed spike. 

We can have models with a clear signature of the spike indicating precisely its 
location because the slope of 5z depends directly on wde and several models are 
displayed in Figure HJ In luminosity-distances the feature of the spike is erased by the 
integration so the effect can easily be degenerate with variations of other cosmological 
parameters and there is no signature which can be unambiguously attributed to the 
presence of a spike in the equation of state parameter wde- The red curve with 
a pronounced spike at z — 0.2 (denoted Spike R in Table 1) is ruled out by SNIa 
data while the green curve with a less pronounced spike located at the same redshift 
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Figure 3: a) The equation of state parameters of our models with spike appearing in 
Figures 2 (same colour as in Figure 2) are shown on the left panel, b) On the right 
panel, the oscillating model used in Figure 2 is displayed. 
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Figure 4: a) SNIa distance modulus data from the Constitution set (blue points) 
and Union2 set (red points), b) Number of supernovae per redshift intervals of 0.1 
for the two data sets. There is a substantial addition in the Union2 data set in the 
interval 0.1 < z < 0.4, which accounts for the stronger constraints in models with 
modifications at low redshift. 



is in tension with the data. On the other hand the model with a spike at z — 1 
(denoted Spike P in Table 1) provides a very good fit to the data. It is interesting 
to remark that, due to the reduced number of data at high redshifts, even a stronger 
modification in Wde like this one can go unnoticed if only SNIa and BAO data is used. 
From Figure 2, the difference in the quantity 10 10 5z between the SpikeP model and 
ACDM at z ~ 2 in a 15-years experiment is estimated to be about 1.2. As things 
stand now this is still too small to allow for a discrimination by a planned experiment 
like CODEX [T4]. However an improvement of order 4 of the sensitivity would allow 
to distinguish both models. 



6 





modi 


mod2 


SSS 


Spike R 


Spike G 


Spike P 


Osc. 


LCDM 




0.267 


0.255 


0.257 


0.241 


0.254 


0.272 


0.262 


0.266 


X2 


544.5 


549.5 


547.0 


554.6 


546.7 


545.0 


544.9 


544.9 



Table 1: Best-fit Q m ,o and x 2 values after marginalization of H , using the Union2 
set [UJ of 557 SN la, the BAO distance ratio [12] and the CMB shift parameter [13] . 
The first three models correspond resp. to the green, pink and red curves of Figure 
1. The models Spike R, Spike G and Spike P correspond resp. to the red, green 
and pink curves of Figure 2a. There are now enough SNIa data on low redshifts 
to put most of the models with a large variation of wde, or a pronounced feature, 
on very low redshifts either in tension with the data (Spike G, SSS) or even to rule 
out some of them (Spike R, mod2). Redshift drift on these low redshifts would be 
virtually impossible to measure. But other models are viable and could be potentially 
constrained using the redshift drift. In particular, Spike P is a strong modification at 
large redshifts while fitting the data almost as well as ACDM. 

3.3 Oscillating dark energy 

Another possibility is provided by models where DE has an oscillating EoS parameter 
(see e.g. [15]), viz. 

Wde = -1 + A osc sin [Sln(l + z) + C] . (16) 

Oscillations in H(z) are induced but will be essentially erased in the luminosity- 
distances. In contrast these oscillations appear clearly in the redshift drift as models 
displayed on Figure [2] show. It is quite clear that a detection of these oscillations 
require very accurate data. These oscillations could be detected with accurate data 
on redshifts z ~ 2 in the range where one expects at least in principle accurate data 
to be possible. 

4 Some models beyond General Relativity 

In this section we would like to show that redshift drift measurements could help in 
some cases to establish whether a model is inside General Relativity (GR) or outside 
it. We have in mind models outside GR for which the gravitational constant defining 
the Chandrasekhar mass evolves with time so that the intrinsic luminosity of SNIa 
would not be constant for such models. We caution that there are many modified 
gravity models where this is not the case because the gravitational constant around 
compact objects is constant and basically equal to the usual Newton's constant (see 
[12] for a related discussion in the context of some viable f(R) DE models). 

For our discussion it is sufficient to consider spatially flat universes. Let us assume 
that the intrinsic luminosity of SNIa has some redshift dependence 

L(z) = L Q(z) Q(z = 0) = 1. (17) 
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In principle this could also arise for physical reasons unrelated to modifications of 
gravity (see e.g.[17j), but we will have in mind DE models where the variation of L 
is induced by the variation of the gravitational constant defining the Chandrasekhar 
mass. For example (II 7p would occur in some models with varying gravitational cou- 
pling (see e.g. [T5]) and in scalar-tensor models. 

From (I17p . the expression for the measured flux T of SNIa can be written as 

-p = L i z ) L ° an 
And 2 L AttV 2 l ' { ' 

where we have introduced the quantity 

V L = d L Q~\ . (19) 



Measuring the flux J 7 does not allow to disentangle the expansion of the universe 
from the variation of the intrinsic luminosity. One can only recover from SNIa data 
the quantity V L (z). 

Confusing the two quantities di and T>i would lead to an erroneous retrieval of 
the background expansion because H~ l (z) ^ y -j+^ J ■ We have instead the equality 



V, 



l + z 



l + z 



(20) 



Measuring H(z), and hence indirectly also di{z), from the redshift drift allows us 
to check whether one has Q ^ 1 and also provides us with information about the 
behaviour of Q{z). 

In the rest of this Section we consider the specific case of (massless) scalar-tensor 
models. In this model, the usual Newton's constant in the equation for the growth of 
perturbations (and in the Poisson equation) is replaced by the effective gravitational 
coupling constant G e s(z), viz. [19] 



h o m +\ — TTz) m = lr ' G m '° dm - W 

Knowing H(z) and the perturbations S m (z), one can check whether these two func- 
tions are consistent within a given model. The quantity G c ff,o = G e s(z = 0) is the 
coupling constant measured in a Cavendish type experiment, its numerical value is 
therefore extremely close to G. 

Taking into account the dependence of the Chandrasekhar mass on G e g and as- 
suming simple SNIa models, the peak luminosity turns out to be proportional to 
' G ^)"i and we have |20] 



Q{z)={^-\\ (22) 
leading to the well-known modification of the distance modulus \x 

fi = 5\ogV L + /i = 51ogrf L + ^log ( —55- J + fi , (23) 

4 V^cff.O/ 
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with Hq = 25 + 5 log 



2Eol 

Mpc 



If the second term is unknown, we cannot retrieve 



H(z) from fx. On the other hand if it is ignored, i.e. identifying T>l with cIl, an 
incorrect H(z) is obtained from the observed \x. While ignoring this modification of fi 
may be enough to show at least an inconsistency with the growth of perturbations in 
ACDM, or more generally in GR with smooth non-interacting DE, it does not allow to 
infer the correct behaviour of the quantity G e s(z) from the growth of perturbations. 
On the other hand, redshift drift measurements yield directly the quantity H(z) 
independently of SNIa data. Combining with SNIa data on a range of redshifts where 
both data overlap, one can recover the behaviour of G e g(z) in this interval applying 
( T2U|) . The growth of matter perturbations give a consistency check with 



5' 2 



\ /2 
% 



m,0 



G e s 5' 5 
G l + zf 



dz' 



(24) 



Additional properties of the growth of matter perturbations can be used as additional 
tests (see e.g. |21j). Finally, let us mention that a varying G e g will also affect the 
width of the light curve which could also provide some information about the evolution 
of G eS . 



5 Summary and conclusions 

We have studied models for which accurate redshift drift data can exhibit distinctive 
features not present in luminosity-distances. We emphasize that another incentive 
to achieve accurate redshift drift measurements comes from some DE models outside 
GR where an evolving gravitational constant affects the intrinsic luminosity of SNIa. 
It is well-known that these models can be efficiently probed by checking whether 
the expansion rate H(z) and the perturbations 8 m (z) are consistent with each other 
inside a given model. The redshift-drift gives a direct probe of H(z) independently 
of the modifications of gravity in contrast to SNIa data and we have illustrated this 
with scalar-tensor DE models. The expansion rate on very low redshifts can also be 
probed with Baryon Acoustic Oscillations (BAO) (see e.g. [22J) which rely on the 
nature and the evolution of matter perturbations. Systematics as well as the probed 
range of redshifts are different, so the two methods would be complementary. It is 
interesting that models with a strong variation of their EoS both on redshifts z < 0.1 
and z ~ 1 would probably escape all observations which leaves some uncertainty 
about the present and future evolution of our universe (see for example ourmodl and 
Spike P). Modifications at high redshifts, in particular, are poorly constrained by 
current data but can be potentially probed by redshift drift measurements. 

We caution that the effect is tiny and that it is not clear yet whether it could 
be measured with the required accuracy. This obviously represents a technological 
challenge. However, as the exquisite data of the Cosmic Microwave Background 
anisotropy have spectacularly shown (see e.g. [13]); technological progress could make 
such measurements possible. If this is the case, redshift drift data could become useful 
as a complementary probe in order to help unveil the nature of DE especially if DE 
is of the kind investigated here. 
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